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Abstract 

For a wide range of models, dark matter can interact with QCD gluons via chromo- 
Rayleigh interactions. We point out that the Large Hadron Collider (LHC), as a gluon ma¬ 
chine, provides a superb probe of such interactions. In this paper, we introduce simplified 
models to UV-complete two effective dark matter chromo-Rayleigh interactions and iden¬ 
tify the corresponding collider signatures, including four jets or a pair of di-jet resonances 
plus missing transverse energy. After performing collider studies for both the 8 TeV and 
14 TeV LHC, we find that the LHC can be more sensitive to dark matter chromo-Rayleigh 
interactions than direct detection experiments and thus provides the best opportunity for 
future discovery of this class of models. 


1 Introduction 


Although there is no doubt that dark matter interacts gravitationally both among itself and with the 
Standard Model (SM) particles, we still have no convincing evidence for other dark matter interactions. 
Among the three forces in the SM, dark matter can not have order-one couplings under the electroweak 
forces. Otherwise, it will either emit photons or scatter off nuclei with a too-large cross section in 
direct detection experiments. Suppressed couplings to the Z boson or the Higgs boson will be probed 
further in upcoming dark matter experiments [T]. For QCD interactions, the dark matter particle 
can not have a color charge due to confinement. If, instead, it is a QCD composite particle, the 
hadronic bound states should have a large interaction strength with SM pions and thus have a large 
scattering cross section with nucleons. This type of dark matter can not penetrate the Earth to reach 
detectors in underground direct detection experiments, but satellite-based X-ray quantum calorimetry 
experiments impose stringent bounds on this scenario [2]. 

Another way for dark matter to interact with gluons is through effective contact interactions, 
comparable to photon-molecule Rayleigh interactions. Introducing such higher-dimensional effective 
operators to describe dark matter chromo-Rayleigh interactions is then independent of whether the 
dark matter is elementary or composite. For instance, in supersymmetric models the neutralino can 
couple to two gluons via a stop-top quark loop. In extra-dimensional models, a similar loop from 
the top quark and its Kaluza-Klein mode can generate the effective chromo-Rayleigh dark matter 
interaction [3]. Composite dark matter models can also generate chromo-Rayleigh interactions, with 
strengths typically suppressed by the composite scale [4]. In this paper, we perform a generic study on 
dark matter chromo-Rayleigh interactions. We pay particular attention to collider searches and the 
potential signatures associated with chromo-Rayleigh interactions. This study is similar to searches 
for electrical Rayleigh interactions of dark matter in Ref. where the electrically charged particles 
can be searched for directly at the LHC. 

For simplicity, we choose the dark matter particle to be a complex scaler denoted as X. We study 
two types of dimension-six interactions: X and (^XX — X^X^'j GG, where G is the QCD field 

tensor. For direct detection searches, the first operator provides spin-independent scattering while the 
second operator provides spin-dependent and momentum-suppressed scattering. For collider searches, 
the standard model independent signature is a mono-jet plus missing transverse energy. Taking into 
account the well-known limitations of such searches [71l24|. we introduce simplified models to UV- 
complete these two operators. Collider signatures owing to such UV-completions depend strongly 
on how the colored mediators decay. We discuss several possible signatures, including multi-jets and 
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pair-produced di-jet resonances with missing transverse energy. 

The rest of this paper is organized as follows. In Section [2l we work out the dark matter thermal 
relic abundance, direct detection cross sections, and collider constraints from the effective chromo- 
Rayleigh contact interactions. In Section [3l we introduce two simplified models to UV-complete the 
previous contact operators. Section [3T] introduces a scalar color-octet with the signature of four jets 
plus missing transverse energy. Sections E^HSS] introduce two fermion color-triplets with the additional 
signature of paired di-jet resonances plus missing transverse energy, and Section 0] offers a comparison 
to direct detection experiments. We conclude in Section [5l 


2 Contact Interactions 


Our model consists of a complex scalar dark matter field, X, which is a singlet under the SM gauge 
group. We introduce the following two CP-conserving, dimension six operators coupling dark matter 
to the gluon field 


^cRayleigh 

^cRayleigh 


dvr 




iOLs 

dvr A 2 


{XX - X^X^)Gl^G^^'', 


( 1 ) 

(2) 


where A* is the cutoff scale and ^ is the dual gluon field strength tensor. The 

overall operator normalization accounts for a loop factor. In order for X to be stable, we impose a 
Z 2 symmetry under which X is odd. Based on these two effective operators, we first calculate the 
thermal relic abundance, direct detection cross sections, and collider constraints. 


2.1 Thermal Relic Abundance 


Depending on the UV physics, the dark matter sector could be more complicated than just one state. 
Therefore, the dark matter thermal relic abundance calculation based entirely on the operators in 
Eqs. m and (I5|) can only provide a guidance for the potential parameter space in Mx and A* for 
thermal dark matter. For the first operator, we have the dark matter self-annihilation rate from the 
process A^A —>• GG as 


{av ){X^X GG) 


1 1 


J 2 

vr^ A| 


= s. 


(3) 


to leading order in the dark matter’s relative velocity v expansion. Here, the overall factor of 1/2 
is due to the relic density being comprised of particles and antiparticles. For the second operator. 
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we have the annihilation rate from the process —)■ GG, in terms of components with 

X = {Xr + iXi)/y/2, as 


- [{av){XRXi 


GG)] = - 




TT' 


3 Xk 


= s. 


(4) 


The dark matter relic abundance is inversely proportional to the annihilation rate and has a formula 
~ 1-07 X 10® GeV“^ Mp\s) with g* as the number of relativistic degrees of freedom at 

the freeze-out temperature and is taken to be 86.25 and the Planck scale is Mpi = 1.22 x 10^® GeV. The 
freeze-out temperature xp is given by Xi? = In [0.05 g MpiMxs/{y/g* xp)] with g = 2 and is typically 
0 ( 20 ). 


2.2 Direct Detection 


For the first operator use the matrix element of inside a nucleon to derive 

the dark matter coupling to two nucleons. The trace anomaly of the QCD energy-momentum tensor 
implies [2511^ 


mN{N\N) = {N\ E + j) -I- 

l< 2 <n f 



asG^G'^nm , 


(5) 


where /3'^f = —(11 — 2nf is the beta function at leading order, Uf is the number of quarks, 

and 7 is the anomalous dimension of the quark field. So, at leading order in as and keeping only 
Uf = 3 light quarks, we have 


{N\asG%G^^‘'\N) = ^mx 


mx 


-{N\muUU + ruddd + msSs\N) — 1 


= Y Kfu + fd + fs) - 1] • (6) 


Recent Lattice QGD updates of the calculation for the strange quark matrix element has fu + fd + fs = 
0.085^0014 [2Z]- In our numerical calculation, we will use (A^|q!sG“j^G“^'^|A^) —2.56mAr. The 

formula of the scattering cross section is then calculated to be 


^si _ 
— 


K m 


N 


dvr Ai {mx + MxY ' 


(7) 


with K.= l{fu + fd + fs-i) ^ -0.20. 

For the second operator, need to know the matrix element of GG inside a nucleon. 

This matrix element is related to the anomalous divergence of the iso-singlet axial current operator 
by 
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In the large-A^c and chiral limit and using the relation {N\ui'^^u+di^^d+si^^s\N) = 0, one finds [28H30] 
{p\^G%G^np) = 389MeV = {n\^G%G^\n) = -2MeV = , (9) 

where the dimensionless parameters are rjp ~ 0.41 and rjn ~ —0.0021 (the instanton calculation in 
Ref. [3T] obtained dramatically different numbers). The large difference between r]p and rjn indicates a 
large isospin violation for pseudo-scalar coupling to nucleons. One can then use the matrix elements 
to translate O 2 to the interaction between dark matter and nucleons, giving 

2z rjNjnN Nij^N . (10) 

In the non-relativistic limit, one has Nij^N ~ 2iq ■ s with q as the exchange momentum of the 
scattering process and s as the spin of a nucleon. So, for this interaction, we have both spin-dependent 
and momentum-suppressed scattering. The spin-dependent differential scattering cross section is 

_ VN 

dcosO 2 7rA| {mj\f + Mx^ 

For Mx ~ 100 GeV, a small cutoff scale of A 2 ~ 100 GeV, and a typical exchange momentum of 
q ~ fJ^XA V ~ 100 MeV, the spin-dependent scattering cross section of dark matter off a proton is 
10“^ pb. This is far below the current direct detection experimental bound |321l33j . 

2.3 Collider Constraints 

For the two operators considered here, the universal signature at the LHG is that of a mono-jet plus 
missing transverse energy mM- For fermion dark matter coupling to two gluons, both the CMS |36j 
and ATLAS [37] collaborations have imposed limits on the cutoffs of the effective operators. To 
estimate the constraints on our scalar dark matter case, we use FeynRules |3S] to generate a model 
file for MadGraph [39]. We then use Pythia |40j to shower and hadronize the parton-level events. 
Finally, we use PGS [H] to cluster hadrons into jets and simulate detector effects. 

Following the same analysis procedure from the CMS collaboration in Ref. |36| at 8 TeV and with 
19.7 fb”^ luminosity, we have found that imposing a cut on the missing transverse energy FI™®® > 
500 GeV (or requiring less than 164 signal events) provides the strongest bounds for a wide range of 
dark matter masses. We show the constraints on the cutoffs of the effective operators in Tabled] for 
different dark matter masses. The constraints stay constant for light dark matter below 100 GeV and 
become weaker as one increases the dark matter mass beyond around 100 GeV. One can also see that 
for heavier dark matter beyond 200 GeV, the constraints become much weaker such that the cutoff 
is even below Mx, which indicates a breakdown of the perturbative description of the effective field 
theory. 
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Mx (GeV) 

Ai (GeV) 

As (GeV) 

I 

130 

170 

10 

120 

180 

100 

120 

180 

200 

110 

160 

400 

90 

130 


Table 1: The collider constraints on the cutoff of the effective operators for different dark matter 
masses at 90% CL. The mono-jet analysis with > 500 GeV from the CMS collaboration in 

Ref. [36] has been used. 

3 Simplified UV-completion Models 

As already can be seen from Table [TJ the constraints on the cutoffs of the contact operators from 
the mono-jet searches are not that stringent. For a 100 GeV dark matter particle, the constrained 
cutoff is just comparable to the dark matter mass. This calls for UV-completed models to reduce the 
uncertainties from an effective field theory description. In this section, we consider several classes of 
models to illustrate that collider signatures beyond the mono-jet may provide a more sensitive probe 
of dark matter chromo-Rayleigh interactions. 

One of the simplest ways to UV-complete the two operators in Eqs. is to introduce a 

color-neutral scalar or pseudo-scalar, which can couple to two gluons either from the top quark or 
new heavy colored fermion loops. For the operator ga,n simply introduce a Higgs- 

portal dark matter coupling like XX^HHK With the effective coupling of the Higgs boson with two 
gluons in the SM, one generates the coupling of two dark matter particles to two gluons. For the 
second operator two-Higgs-doublet-portal dark matter models have a pseudo-scalar as 

a mediator to generate the dark matter chromo-Rayleigh interactions (see Refs. [T9l[20l|42l|43] ). The 
existing search strategy for the Higgs-portal or pseudo-scalar-portal dark matter should cover this 
class of UV-completion models jUj. Therefore, we do not discuss this case in great detail. Instead, 
we consider two other ways to UV-complete the two operators and point out more interesting collider 
signatures. 

3.1 QCD-charged Particle Mediation for 

In this subsection, we examine a class of models with additional QCD-charged scalars to generate the 
effective chromo-Rayleigh interaction of gpggifjgj^y’g gake, we introduce a real color- 
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Figure 1: Representative loop diagrams for generating O 


cRayleigh 

1 


via colored scalar loops. 


,0- 


octet and electroweak-singlet scalar, LU with a = 1, • • • ,8 and a mass of which, for example, 

appears in the Renormalizable Coloron Model (ReCoM) |45p46j. At the renormalizable level, one has 
the following interaction coupling two Gj^’s to two dark matter helds 


£ D G\X^X 


( 12 ) 


At one-loop level, the effective operator can be generated through the diagrams in Fig. [T] and has 
the following calculated form 

A a. 


F{t) 


8M2^ 




(13) 


where the form factor F{t), with r = {pi + p 2 )‘^/{4:Mq^), is one in the limit r —)• 0. Here, pi and p 2 
are the gluon momenta pointing to the vertex. In the limit that Mq^ S> Mx and with small gluon 
momenta much below \pi +P 2 \, one can match the colored particle mass to the effective operator cutoff 
defined in Eq. ([T]) to obtain 


a2 _ 8 2 

^'■1 - ■ 


(14) 


For this class of UV-complete models, one can also work out the collider constraints. Depending 
on whether the color-octet is odd or even under the dark matter Z 2 , one has different decay channels 
for G^. If it is Z 2 -eveii, the following renormalizable Gh cubic self-interaction can make Gh unstable 

/^G dabc GhG^h^h I (15) 

with dabc the totally-symmetric 5'C/(3 )qcd tensor. This introduces the decay of Gh into two gluons 
at one loop with a width of [45|I17] 

15 


F(G^^ ^ gg) = 




9 


( 16 ) 


similar analysis can be performed for other QCD representations. 
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-Z 2 -even Gh can be pair-produced at the LHC from their QCD interactions, which then decay into 
paired dijets. Following the same reinterpretation of the experimental data as in Ref. [48], we find 
that the current searches for paired dijet resonances at the 8 TeV LHC have set a constraint on its 
mass of [l9| 

Mgh > 520 GeV, for Za-even Gh , (17) 


where we have only included the QCD productions. 

If Gh is Z 2 -odd, the dark matter particle X has to appear in the Gh decay products. The 
operators mediating Gh decaying to X hrst occur at the dimension six level and contain the following 
two parity-conserving operators 


-h h.c. 


G%XHQLtHR 


+ h.c.. 


(18) 


A'2 ' ■ A" 2 

Here, with a = 1,2,--- ,8 are the 5'17(3)qcd generators. For the second operator, we have only 
included the top quark by assuming that the coupling is proportional to the quark mass. Using the 
equation of motion, = —Qs the first decay-inducing operator becomes 


^a%d,xi}jVq. ( 19 ) 

As a result, the decay channels of the colored state Gh are mainly Gh Xqq. For each flavor, the 
three-body decay width is calculated to be 


r{GH 


9 s 


Xqq) = 

q{MG„-Mx? 


3 • 29 tt^A'i 


MlVlc 

1 J 



Mi 


+ (M2^-s)2-2M2(M2^+.) 


m2 


3/2 

( 20 ) 


Choosing gg = 1.1, = 500 GeV, Mx = 10 GeV, A'^ = 1 TeV, and neglecting the quark mass, 

F = 0.0002 GeV for each flavor and Gh can decay promptly for collider studies. 

From the above assumptions, the collider signature for pair-produced Z 2 -odd Gh contains 4j -|- 
and 2t-|-2t-|-Fi™®®. Generically, the top-quark rich final state can be easily searched 
for at the LHC. So, we concentrate on the first operator in Eq. (1181) and derivate a more conservative 
bound based mainly on the 4j -|- E™'®® hnal state. Following the analysis in Ref. [SD] at the 8 TeV 
LHC with 19.5 fb“^, we find that the set of cuts with 300 < fix < 450 GeV, 800 < Hx < 1000 GeV, 
3 < Vjets < 5 with px{j) > 50 GeV and |r/| < 2.5 provide the best constraint. For Mg^ = 600 GeV 
and Mx = 10 GeV, we have Br(GH —)■ 2j + X) = 97.5% and Bt^Gh tt + X) = 2.5%. The 
production cross section at the 8 TeV LHG after this set of cuts is approximately 1.9 fb x K with 
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the ii'-factor as 1.8 m- This amounts to a total of 65.4 events at 19.5 fb which is very close to 
the allowed number of signal events (65.6 from Ref. [50]) at 90% C.L. So, for a large mass splitting 
between the color-octet state Gh and the dark matter X, the current LHC bound is 

Mgh ^ 600 GeV, for Zs-odd Gh ■ (21) 

Independent on how Gh decays, we have set a constraint on the Gh mass to be above 500-600 
GeV. Using the relation in Eq. (1141) . this constraint translates to a limit on the effective cutoff scale 
of Ai > 1.5 — 1.7 TeV for A = 1, which is far more stringent than the limits in Table [1] set from 
the mono-jet search. To compare with the limits from direct detection experiments, we show the 
reinterpreted collider limits using Ai > 1.7 TeV for the first chromo-Rayleigh interaction in the left 
panel of Fig. [7| The left panel of Fig. [7] shows a comparison between the limits on the scattering cross 
section obtained from direct detection, monojet, and Aj + F™'®® searches. 


3.2 QCD-charged Particle Mediation for 

To UV-complete the operator introduce the following CP-conserving Lagrangian 

£ D -yi (^X -b (V’iV ’2 + {y 2 '0i75V’2 + y 2 V'275V’i) • (22) 

Here, ipi and 'tp 2 are chosen to be QCD triplets and their electroweak quantum numbers will be 
discussed and specified later. To conserve the dark Z 2 symmetry, one of the fermion triplets must 
be F 2 -odd while the other is F 2 -even. Using the freedom of field redefinitions of tpi and ^ 2 , one can 
keep the first coupling, yi, to be a real number while the second coupling, 2/2 j aiay in general be 
complex. We choose 2/2 to be a real number to satisfy GP symmetry. In terms of the components, 
X = {XR + iXi)/V2, we have the interactions 

C D -V 2 yi Xr {i)iil)2 + i’ 2 '^ 1 ) - V2iXi {y2 '0i75V’2 + 2/2 V'275V'i) ■ (23) 


To preserve G and P, one has Xr to be C-even and P-even, and Xj to be C-even and P-odd. In 
terms of Xr and X/, we have the effective operator 


O 


cRayleigh _ 


a.9 


27rA2 


\r \r /^a /^a 

JL I ij 




(24) 


which conserves both C and P. 

At one-loop level, one has the box diagrams in Fig.[2]to generate the effective operator 
the heavy particle limit with ^ Mxj^iMxj , we have the matching condition (see Appendix lAl 

for a more detailed calculation) 

' “ I2/12/2I 
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(25) 




ysmoQSu 


XR{ki) G“^(pi) ^SiQQQQSU 


XR{k,) 


g''"(p 2 ) vQQQQQW 


^2(g) 
i I V'2 V'G ^ 
i’2 


Xi{k2) G^-{P2) <smQQSO 


V'2(g) 

i L '02 01 ' ' 

02 


' X,(fc2) 


Figure 2: Representative loop diagrams for generating yjg^ colored fermion loop. 


For this specific UV-completion model, the first chromo-Rayleigh interaction can also be generated, 
since it does not break any discrete symmetries. 

If Im(yiy 2 ) 7^ 0, the Lagrangian in Eq. (I23p is P-conserving but C-breaking, so CP is also broken.]^ 
The lowest dimensional effective operator that is C-odd and P-even is at dimension-10, for instance 
XrXj , further suppressed by powers of the cutoff scale. If parity was 

broken, lower-dimensional operators like XrXiG’^^G°‘^'^ could be generated with a stringent bound 
from the neutron electric dipole moment. 

As mentioned before, in order to conserve the dark Z 2 symmetry the two fields, 0i and 02, should 
have opposite dark parity. Fixing the mass relation we have two cases. The case A 

has Z 2 -even 0i and Z 2 -odd 02 and the case B has Z 2 -even 02 and Z 2 -odd 0i. For both cases, we 
introduce the following dimension-5 operator to mediate the 2^2-even particle decaying into two jets, 

(26) 

Here, i = 1(2) for case A(B); the electroweak quantum numbers of 0j are then the same as ur (or 
dpi). As a result of this decay mode, the unstable particle 0i behaves as a dijet resonance at colliders. 
One could also consider other potential decay channels by introducing the electroweak dipole moment 
operator, which has a clearer signature at colliders. 

For case A with a 2^2-even 0i, the collider signature of 0i is pair-produced dijet resonances. The 
searches at the 8 TeV LHC have set a constraint on its mass [l9] as > 500 GeV. Although the 
signature is interesting by itself, the discovery of this dijet resonance can not prove that dark matter has 
been produced at colliders. For the Z 2 -odd 02 particle, ithe decay channel is 02 —>■ X -|-0i — >• A -|- 2j. 
Depending on the mass splitting of 0i and 02, one has 0i to be off-shell for < Mx and 

on-shell for > Mx- For the off-shell intermediate 0i case, this signature is very similar 

to the SUSY squark searches with a heavy gluino [50], except with a larger production cross section 

^For a complex j /2 with CP-violating interactions, one may wonder about generating the Weinberg operator |52 |. 
fabcG’^pG^PCk^^''. We note that the P conservation forbids the generation of the effective Weinberg operator. 
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than a single flavor squark. For dark matter mass Mx ^ 100 GeV, the constraint on the ^^2 mass is 
^-^2 ^ ^50 GeV, assuming that the signal acceptance is similar to the squark one [50]. For a heavy 
dark matter mass close to 290 GeV, the constraint becomes weaker and is ^ ^00 GeV. Using 
the conversion formula in Eq. (1251) , we show the constraints from 4j + on the effective cutoff for 
different dark matter masses in the left panel of Fig. [S] by setting — Mx and yiy 2 = 1- 




Figure 3: Left panel: constraints on the cutoff of the effective operator fgp Q^se A with 

Z 2 -even 'ipi and Z 2 -odd V’ 2 - Right panel: constraints for case B with Z 2 -odd 'tpi and 2)2-even '02- 


For the on-shell 'ipi case, the signature is more interesting with a pair of dijet resonances plus 
a large The relevant production Feynman diagram is shown in the left panel of Fig. [4| We 

perform a detailed collider study for this interesting signature in Section 13.31 As a comparison, in the 
left panel of Fig. [3] we show the interpreted constrains for Mx ^ 100 GeV on the matched cutoff from 
Eq. (|25p for both the 8 TeV LHC with 19.5 fb“^ and the 14 TeV LHC with 100 fb“^. One can see a 
clear improvement on constraining this effective chromo-Rayleigh interaction from a dedicated search 
beyond the simple mono-jet signature. 

For case B with Z 2 -odd ijji and Z 2 -even "^ 2 ; the searches for pair-produced dijet resonances provide 
a constraint of > 500 GeV. The decay of ipi needs to go through an off-shell ^^2 via ipi ^ X+11^2 
X + 2j. Applying the SUSY squark searches with a heavy gluino |50|, the constraint on the i/ji mass 
is > 850 GeV for Mx ^ 100 GeV and > 500 GeV for Mx close to 290 GeV. In the 

right panel of Fig. [H] we show the reinterpreted constraints on the effective cutoff for yiy 2 = 1. 
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3.3 Pair-produced Dijet Resonances plus 

In this section, we perform a detailed collider study for the process of pp V’ 2 V ’2 X+il^i+X "^+^ 

+ 4j. The collider signature is pair-produced dijet resonances plus missing transverse energy. In 
the left panel of Fig. [4l we show the Feynman diagram for this process. The QCD production cross 



xt 



Figure 4; Left panel: the representive Feynman diagram for the process of pp —)• ip2'4’2 X ipi -\- 
Xl + -4)1^ XXl -|- 4j. Right panel: the tree-level production cross sections of pp 'ip 2 ' 4 ’ 2 - 


sections for 'i/’ 2 V ’2 at the LHC are the same as vector-like t't' 153]. We show the tree-level production 
cross sections, calculated using MadGraph [39], in the right panel of Fig. [4| In the following analysis, 
we will ignore the signal X-factor because we will use tree-level cross sections for backgrounds. □ 

The main SM background comes from ZjW^ + n jets with Z ^ vv and leptonic decays of W^. 
After comparing the results from matching parton showers and matrix elements [54] . we have found 
that the background of parton-level ZjW^ + 4 jets with ptU) > 120 GeV provides a good estimation 
of total ZjW^ + n jets background for the 8 TeV LHC. Therefore, we use ZjW^ + 4 jets as an 
approximation to save simulation time. There also exist additional, sub-dominant semi-leptonic ti 
backgrounds, which will be kept in our analysis. 

Starting from the 8 TeV LHC with 19.5 fb“^, we choose the model point {m ^-^, = (600,900) GeV 

as a benchmark to optimize our cuts on kinematic variables. We choose the basic cuts on the jet and 
missing transverse momenta to be ptUz) > 140 GeV and > 275 GeV and required at least four 
jets satisfying the jet px cut in the final state. To reduce the tt and VF^-t-jets backgrounds, we also 
veto events containing a lepton with pt{(-) > 20 GeV. Since our signal has a pair of dijet resonances 
with the same mass, we choose the combination among three possible dijet pairs with the smallest 

^One might expect the overall significance in the later analysis to be increased by a factor of ~ 1.2 — 1.3 due to the 
inclusion of NLO effects. 
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?/^avg (GcV) 


Figure 5: Left panel: the distributions of the fraction of events as a function of Arndyet/^avg 
for the signal of = (600,900) GeV and the summed background. Right panel: the 

averaged dijet invariant mass distributions for the signal and various backgrounds, after the cut of 
Amdijet/m-avg < 0.4. 


dijet invariant mass difference, Amaijet- To further reduce the SM background, we show the event 
fraction histogram distribution in terms of the variable Amaijet/^avg in the left panel of Fig. [5l Signal 
events prefer a smaller value of dijet invariant mass difference than background events. As a result, 
we impose a cut on this variable with Amaijet/nravg < 0.4 to further increase the discovery sensitivity. 
In the right panel of Fig. [U we show the averaged dijet invariant mass distribution for the signal and 
backgrounds. Since the signal events mainly distribute around the 'ipi particle mass, we also impose 
one additional cut with |mavg — < 0.1 to further improve the sensitivity. With the above 

cuts, we show the values of S/y/B for different mass points in the left panel of Fig. [H At 90% C.L., our 
simulated results show that the model point of = (650,950) GeV can be covered, which 

corresponds to A 2 > 1.1 TeV. Once again we see a dramatic improvement relative to the standard 
mono-jet signature. 

At the 14 TeV LHG with 100 fb“^, we choose two benchmark points with Mx = 10 GeV and 
= (900,1800) [(900,1500)] GeV. To optimize the discovery sensitivity, we choose the 
following set of cuts: at least four jets with pxiJ) > 200 GeV; > 900 [600] GeV; Amdijet/w-avg < 
0.4 [0.3]; Im-avg —wii/'il < 0.1m.i/,j. The signal significances, S/y/B, for different model points are shown 
in the right panel of Fig. [6l taking the larger S/y/B value from between the two benchmark points. 
We show the combined 90% G.L. exclusion contour in the red dashed line, which has the lower right 
region covered by the set of cuts for (900,1800) GeV and the upper right region by (900,1500) GeV. 
At 90% C.L., one can see that the model point of ( 771 ^^, 771 ^ 2 ) = (950,1800) GeV can be excluded, 
which corresponds to A 2 > 1.8 TeV. 
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Figure 6: Left panel: the signal significance S/y/B for different model points with the cuts optimized 
for {m^^, 111 ^ 2 ) — (600, 900) GeV and Mx = 10 GeV at the 8 TeV LHG with 19.5 fb“^. The red dashed 
line is the extrapolated 90% G.L. exclusion curve. The black dashed line has • R-ight 

panel: the same as the left panel but for 14 TeV with 100 fb“^ and the set of cuts for the combined 
benchmark points { 171 ^^, 171 ^ 2 ) — (900,1800) GeV and (900,1500) GeV with Mx = 10 GeV. 

4 Comparison to the Direct Detection Limits 

To compare with the constraints from direct detection experiments, we show the interpreted dark 
matter-nucleon scattering cross sections from the LHG mono-jet, 4j -|- S™'®® and paired dijet-|-£1™‘®® 
searches in Fig. [71 In the left panel, we show the spin-independent scattering cross section including the 
constraints from the LUX collaboration [36]. Gompared to the direct detection limits, the constraints 
from the LHG are not strong for heavier dark matter masses, but are more stringent for a light dark 
matter mass below around 5 GeV. Compared to the limits from mono-jet searches, the 4j -|- £^™'®® 
signature is definitely a more sensitive channel for constraining the spin-independent scattering cross 
section from the first chromo-Rayleigh interaction, 

In the right panel of Fig. [71 we show the interpreted collider constrains for the second operator, 
Here, we only show the dark matter-proton scattering cross section, since it 
has a larger value than dark matter-neutron one for the same dark matter mass. For the interpreted 
limits from collider searches and because of the additional momentum suppression with q ~ fixAV, we 
fix the dark matter velocity to be u = 10“^ and choose a typical target nucleus mass mx = 100 GeV 
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in Eq. m- Compared to the constraints from direct detection experiments, the collider searches 
obviously provide a more sensitive probe of this type of dark matter chromo-Rayleigh interactions. 




Figure 7: Left panel: the constraints on the dark matter spin-independent scattering cross sections 
from the LHC mono-jet [36] and multi-jet |50| searches at the 8 TeV LHC with ~ 20 fb“^ 
and the LUX direct detection experiment |55|. Right panel: the same as the left one but for spin- 
dependent dark matter-proton scattering cross sections from the LHC, SIMPLE [56], PICASSO |57j . 
PICO-2L [58], IceCube [59| and Super-K |^. The scattering cross section is suppressed by both the 
spin-dependent scattering and the exchanging momenta. 


5 Discussion and Conclusions 

One advantage of introducing simplified models is that the relevant model is complete from the renor- 
malizability point of view. There is no worry about trusting theoretical descriptions of the collider 
limits like in the effective operator approach. On the other hand, the simplified models typically 
contain more model parameters and make the comparison with the results from dark matter direct or 
indirect detection less model-independent. For the simplihed models studied here, the collider signa¬ 
tures highly depend on how the QCD charged particles decay. For the first operator, the color-octet 
scalar can decay into two gluons or two quarks plus the dark matter particle. For the second operator, 
the two color-triplet fermions could have many other possible decay channels beyond the ones in our 
paper. For instance, one could replace the operator in Eq. (|26]) by , which can lead to the 

decay of ijj u + '^jZ. The final collider signature could be 2j -|- 2'y{Z) + F™'®® with or without a pair 
of j + j(Z) resonances. To cover the majority of simplihed models for dark matter chromo-Rayleigh 
interactions, one should search for a wide range of potential signatures. 

In conclusion, we have studied simplihed models to UV-complete the chromo-Rayleigh interactions 
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of dark matter. For the first operator, a new color-octet scalar particle could be in reach of the LHC. 
It may decay into just two jets or two jets plus missing transverse energy. The 8 TeV LHC can 
already constrain this color-octet scalar mass to be above 500-600 GeV, which can be translated into 
a constraint on the cutoff of the effective operator to be above 1.5-1.7 TeV. For the second operator, 
two QCD-charged fermions are predicted for the simplified model in this paper. A collider study of 
pair-produced dijet resonances plus missing energy at the 14 TeV LHC can constrain the geometric 
mean of the two fermion masses above around 800 GeV, which can be interpreted as a bound on the 
effective operator cutoff above around 1.8 TeV. 


Note added: We note here that during the completion of our paper, another paper m appeared 
containing some overlap with our UV-completion of the hrst operator, . 
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A Loop Level Calculations for the Coefficient of 

Starting from the Lagrangian in Eq. ()23p . we match the coefficient for one-loop level. 

Other than the two Feynman diagrams in Fig. [2l we also have another four diagrams in Fig. [8] and 
Fig. [9l 











Xi{k2) G''"(p 2) vQQQQiKD 



, XRih) 


" Xi{k2) 


Figure 8: Loop diagrams for generating ^ colored particle in loop.. 


Summing the two diagrams in Fig. [2l we have the matrix element as 


iM-i 


= 4iffsyiy2Tr[f“t^]e^(pi)e,,(p2) x / 


(27rF [q2_rn‘j][(^q-pi)'^-m^][{q-ki-k2)'^-ml][{q-ki)'^-rnf] 


• (27) 
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XR{h) 


XRih) 



In the heavy fermion limit, we have a simplified version of the matrix element as 


iM, = z i^5'^'’e^{pi)e,{p2)e^^’^>’yipP2o 


—nil + 2 m\m 2 log + mim^ 


ni2{nii — m^) 


2\3 


(28) 


The two diagrams in Fig. [8] have a similar answer as zA ^2 = -H- m 2 ). The two digrams in 

Fig. [9] has a summed matrix as 




{pi)£:iy {P2 ) e^''^'"pipP2o 


j mim 2 [2(mf — m^) — (m^ + m^) log (mf/m^)] 


1 


{mf — 


(29) 


Add all diagrams together, we have the summed matrix element 

2 

iMi + iM2 + iMs = 5°-^e^{pi)eu{p2)e^''^''pipP2a ■ 

dvr^ mi m 2 

To match to the coefficient of one has 

as yiy2 _ 2 , Us 
vr mi m2 TT A2 


(30) 


(31) 
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